WHAT IS CLAIMED IS: 

1. A method for approximating y(n)=l/x(n) in FM demodulation, 
where x(n)= I 2 (n)+Q 2 (n), comprising: 

(a) receiving a prior estimated value of l/x(n); 

(b) receiving a present value of x(n); 

(c) adjusting the prior estimated value of l/x(n) to compensate 
for an error between the prior estimated value of l/x(n) and the present value 
of l/x(n); and 

(d) outputting the adjusted prior estimated value of l/x(n) as the 
present value of T/x(n). 

2. The method of claim 1, wherein the prior estimated value of 
l/x(n-l) equals l/(I 2 (n-l)+Q 2 (n~l)), wherein I(n) is an input signal and Q(n) is 

-■ 

a quadrature-phase signal of the input signal I(n). 

3. The method of claim 2, wherein the input signal I(n) comprises 
a band pass filtered secondary audio program signal. 

4. The method of claim 1, wherein the present value x(n) equals 
I (n)+Q (n), and wherein I(n) is an input signal and Q(n) is quadrature-phase 
signal of I(n). 

5. The method of claim 4, wherein the input signal I(n) comprises 
a band pass filtered secondary audio program signal. 

6. The method of claim 1, wherein an error signal equals (1- 
x(n)y(n-l))a, wherein x(n) = I 2 (n)+Q 2 (n), y(n-l)= l/(I 2 (n-l)+Q 2 (n-l)) 5 I(n) is 
an input signal, Q(n) is a quadrature-phase signal of the input signal I(n), and 
"a" is a scaling coefficient/ 
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7. The method of claim 6, wherein the input signal I(n) comprises 
a band pass filtered secondary audio program signal. 

8. The method of claim 1, wherein the Y(n) signal equals y(n-l) + 
(l-x(n)(y(n-l))a, wherein x(n) = I 2 (n)+Q 2 (n), y(n-l)- l/(I 2 (n-l)+Q 2 (n-l)), I(n) 
is an input signal, Q(n) is a quadrature-phase signal of the input signal I(n), 
and "a" is a scaling coefficient. 

9. The method of claim 8, wherein the input signal I(n) comprises 
a band pass filtered secondary audio signal. 

10. A method for demodulating an FM signal FM(n) from a 
secondary audio program signal, comprising: 

(a) receiving in-phase I(n) and quadrature-phase Q(n) portions 
of the FM(n) signal 

(b) generating a first portion of the FM(n) signal that is equal to 
I(n)Q(n)-I(n)Q(n); 

(c) determining a value z(n) based on the first portion of the 

FM(n) signal; 

(d) generating a second portion of the FM(n) signal that is 

9 9 9 9 

equal to 1/ I (n)+Q (n), wherein I (n)+Q (n) is equal to x(n) and y(n)=l/x(n); 

(e) generating a value for y(n) based on l/x(n) that equals y(n- 
l) + (l-x(n)y(n-l))a; and 

(f) multiplying the z(n) value and the y(n) value to produce the 

FM(n) signal. 
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11. A system for approximating y(n)=l/x(n) in FM demodulation, 

2 2 

where x(n)= I (n)+Q (n), comprising: 

means for receiving a prior estimated value of l/x(n); 

means for receiving a present value of x(n); 

means for adjusting the prior estimated value of l/x(n) to 
compensate for an error between the prior estimated value of l/x(n) and the 
present value of l/x(n); and 

means for outputting the adjusted prior estimated value of 
l/x(n) as the present value of l/x(n). 

12. The system of claim 11, wherein the prior estimated value of 
l/x(n-l) equals l/(I 2 (n-l)+Q 2 (n-l)), wherein I(n) is an input signal and Q(n) is 
a quadrature-phase signal of the input signal I(n). 

13. The system of claim 12, wherein the input signal I(n) comprises 
a band pass filtered secondary audio program signal. 

14. The system of claim 11, wherein the present value x(n) equals 
I 2 (n)+Q 2 (n), and wherein I(n) is an input signal and Q(n) is quadrature-phase 
signal of I(n). 

15. The system of claim 14, wherein the input signal I(n) comprises 
a band pass filtered secondary audio program signal. 

16. The system of claim 11, wherein an error signal equals (1 - 
x(n)y(n-l))a, wherein x(n) = I 2 (n)+Q 2 (n), y(n-l)= l/(I 2 (n-l)+Q 2 (n-l)), I(n) is 
an input signal, Q(n) is a quadrature-phase signal of the input signal I(n), and 
"a" is a scaling coefficient. 

17. The system of claim 16, wherein the input signal I(n) comprises 
a band pass filtered secondary audio program signal. 
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18. The system of claim 11, wherein the Y(n) signal equals y(n-l) 
+ (l-x(n)(y(n-l))a, wherein x(n) = I 2 (n)+Q 2 (n), y(n-l)= l/(I 2 (n-l)+Q 2 (n-l)), 
I(n) is an input signal, Q(n) is a quadrature-phase signal of the input signal 
I(n), and "a" is a scaling coefficient. 

19. The system of claim 18, wherein the input signal I(n) comprises 
a band pass filtered secondary audio signal. 

20. A method for approximating y(n)-l/x(n) in FM demodulation, 
where x(n)^ I 2 (n)+Q 2 (n), comprising: 

(a) receiving l/x(n-l); 

(b) receiving x(n); 

(c) adjusting l/x(n-l) to compensate for an error between 
l/x(n-l) and l/x(n); and 

(d) outputting the adjusted l/x(n-l) as l/x(n). 
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